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Abstract. Electrostatic interaction between ions in an ionic channel and the charge 
fiuctuations in the channel mouth are considered. It is shown that the charge fluctuations 
can be enhanced in the channels with low dielectric constant and result in strong modulation 
of the potential barrier at the selectivity site. The effect of charge fluctuational on transition 
probabilities in other molecular dynamical systems is briefly discussed. 



1. Introduction 

Electrostatic interactions and fluctuations play crucial role in controlling transition probabilities, 
conduction and selectivity in nanoscale molecular systems. It is also well- know that the long- 
range Coulomb interaction leads to the strong correlation of ions motion and is responsible for 
the ions distribution and the energy fluctuations in electrolytes [HI [7]. However, in traditional 
approach to the calculations of the escape rates in molecular biology the charge fluctuations 
in the bulk are usually neglected [6]. In this paper the effect of charge fluctuations on the 
probabilities of transition of ions through an open ionic channel is considered. 

Specifically, we consider ion transport across lipid membrane with small dielectric constant. 
The motion is described as a self-consistent solution of the coupled Poisson and Langevin 
equations for ions moving in an open ionic channel [HI H] . The distribution of ions arrival time 
and the charge fluctuations at the channel mouth are investigated numerically and compared 
to the theoretical predictions. Next, the electrostatic energy profile of one ion in an aqueous 
channel through a lipid membrane charged at the selectivity site is calculated. Finally, we 
consider the electrostatic coupling between the charge fluctuations in the channel mouth and 
the energy fluctuations for the ion at the selectivity site. 

It is shown that the volume charge fluctuations in the channel mouth can be modeled as a 
generalized short noise and result in the strong modulation of the potential barrier for an ion 
at the selectivity site at the sub-nanosecond time scale. It is further demonstrated that the 
enhanced by the channel modulation of the potential barrier is much larger as compared to the 
energy of Coulomb interaction between ion in the channel mouth and the ion at the selectivity 
site. We therefore argue that open ionic channels can be viewed as electrostatic amplifiers 



of the charge fluctuations. An analytical approximation to the effect of amplification of the 
electrostatic interaction between the ions in the channel is provided. 



2. Method 

The system considered is made of three compartments of equal size. The middle block constitutes 
the protein through which there is a cylindrical hole representing an approximation of an open 
channel. The channel is bathed by a solution of NaCl on its left and right ends. The electrostatic 
force and potential are produced by the moving ion and fixed charge at the surface, particularly 
the net negative charge bared by a fraction of the protein at the channels selectivity filter. The 
Poisson equation is solved in 3D space with cylindrical symmetry. The electrical sources are 
the fixed negative charge on the protein and the ions located on the channel axis. Since the 
dielectric value of the aqueous pore is not yet known for narrow channels, as first approximation, 
we use the same dielectric value as for the aqueous bath e2 = 80. The dielectric value of the 
protein is taken to be equal to ei = 2. 

The motion of the ions is modeled within a self-consistent framework of Brownian dynamics 
(BD) coupled to the Poisson equation. 
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where rrii, Xi and qi are the mass, position and charge of the ith ion. In the Brownian 
dynamics simulations, water molecules are not taken into account explicitly and are represented 
in Eq. ([T]), by an average frictional force with a friction coefficient rriiji and a stochastic 
force ^/2mi^ikBT^i{t) arising from random collisions. The long range Coulomb interaction is 
represented by the 1/r potential. The addition of the pairwise repulsive soft-core interaction 
potential insures that ions of opposite charge, attracted by the inter-ion Coulomb force, do not 
collide and annihilate each other. r,j is the distance between ions i and j. Uq and Rc are 
respectively the overall strength of the potential and the contact distance between ions pairs [8j . 

We used the following parameters for the simulations: 

Dielectric constants: £2 = 80, £1 = 2; 
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Masses (in kg): niNa = 3.8 x 10"^^ mci = 5.9 x 10 
Diffusion coefficients (in m^s~^): D^a = 1-33 x 10^^, Dci = 2.03 x 10"^ 
(Note that D is related to the friction coefficient via D = ^^); 

Ion radii (in A): r^va = 0.95, rci = 1.81; 
Temperature: T = 298 K. 



Ions 


Uo(kBT) 


i?c(A) 


Na-Na 


0.5 


3.50 


Na-Cl 


8.5 


2.76 


Cl-Cl 


1.4 


5.22 



Table 1. Parameters used in the calculation of the short range ion- ion interaction 



The dielectric force acting on the ion as it moves on the channel axis F^h is calculated 
numerically by solving Poisson equation 

- V- (e(f)V0(f)) = i5 = eE, E = -Vcj) (3) 
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using finite volume method (FVM) In Eq. ([3]), e = £^£0 is the space dependent dielectric 
function with A; = 1 or 2. Sq = 8.85 • 10~^'^ CV^^m^^ is the dielectric constant of empty space. 
p is the source density, (j) is the potential, D and E are respectively the displacement vector 
and the electric field. An effective dielectric constant is introduced at the interface between 
water and the protein. The later procedure is more appropriate when the dielectric function 
present severe jump. In order to preserve the system symmetrical symmetry, the moving ion is 
bound to move on the channel axis. Dirichlet boundary condition is used to fix the value of the 
electrostatic potential at the left and right boundaries of the system as used in experimental 
measurement. Neumann boundary condition is used to set the value of the normal component of 
the electric field. Standard iterative method is used to solve algebraic linear system grown from 
discrete Poisson equation. The results are stored in tables allowing for the fast self-consistent 
BD simulations. 

3. Distribution of arrival times and the charge fluctuations 

In this research we are focused on estimating the effect of the charge fluctuations on the transition 
probabilities of the ions through an open ionic channel. To this end we would like to calculate 
how much time ions spend in the channel mouth and how often do the arrive to the channel 
mouth. Estimation of arrival time of ions at the channel mouth from a solution assuming charge 
neutrality and no applied field can be obtained by considering pure diffusion of ions through 
the hemisphere at the mouth of the channel and is given by Tar = l/(27rzjHZ)jrco); where r is 
the channel radius, cq the bulk concentration, H the Avogadro number, Zi and Di the charge 
valence and the diffusion coefficient of the ith ion (BJ ?]. For an NaCl solution of concentration 
Co = 500mM, estimated arrival time Tar ~ 2.9 x 10~9ses for Na"*" and Tar ~ 3.8 x 10~9ses for 

cr. 

To estimate the charge fluctuations from the simulations we have recorded continuously 
(during a few microseconds) the total positive and negative charge in the channel mouth with 
volume vm = T^f"^ * f, where r = QA. The arrival time was estimated by recording the interval 
of times between the subsequent events of ion arrival to the channel mouth. The results of 
simulations are summarized in the Fig. [TJ It can be seen from the figure that the arrival time 
distribution follows exponential distribution with mean arrival time Tar ~ 3.6 x 10~9ses for Na"*" 
and Tar ~ 4.7 X 10~9ses for Cl^. 

This corresponds corresponding to a Poisson process as expected from the theory. The 
corresponding charge fluctuations at the channel mouth is represented on Fig. [TKb). 

4. Amplification of the electrostatic interaction between two ions in the channel 

We now consider the effect of the charge fluctuations on the transition probabilities of the ions 
through an open ionic channel. First, we analyze the energy profile faced by one ion moving 
through the channel. 

The total electrostatic energy required to move one ion from the bulk solution to a bare 
channel that is 30 Along and 4 Ain radius as a function of the position of the ion is calculated 
as shown in the Fig. [2](a). The potential drop across the channel is A<I> = 90 mV, the charge 
at the selectivity site (at the middle of the channel) is -le. Note that the value of the potential 
barrier for the ion at the selectivity site to exit the channel is AEq as shown in the Fig. [2] (b). 
We now consider the reduction of this potential barrier induced by the second ion situated at 
the channel mouth. The total energy of the channel as a function of the position of the first 
ion moving along the channel when the second ion is located at the channel mouth is shown in 
the Fig. EKb) by dashed line. It can be seen that the presence of a second ion at the channel 
entrance decreases the energy barrier to AEi. We note the reduction of the potential barrier 
from AEq to AEi is much stronger in the channel with low dielectric constant. I.e. the effect 
of charge fluctuation is strongly enhanced by the protein low dielectric constant compared to 
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Figure 1. (a) The arrival time distributions for positive (red line) and negative (blue line) ions 
for cylindrical channel of radius r = 4 A. (b) Fluctuations of the charge at the mouth entrance. 




Figure 2. (a) The potential energy profiles in a cylindrical channel of radius r = 4 A when an 
electric field of 10^ V/m is applied in the z direction. The solid (blue) and dashed (red) lines 
correspond to the channels with and without fixed charges, respectively. The profile of a passive 
channel {eprotein = 80) is indicated by the dotted (black) line. The channel is situated between 
30 and 60 A. (b) Total energy of the channel as a function of the position of the ion when: 
the first ion is located at the channel mouth (dashed line); the channel is empty (solid line); no 
channel (dotted line). Vertical dashed-dotted lines show channel entrance. The height of the 
potential barrier for the ions at the selectivity site and no ion at the channel mouth is denoted 
AEq. In the presence of the second ion in the channel mouth this barrier is reduced to AEi. 



water, in the absence of the protein walls the interaction between two ions will be purely due 



to the Coulomb forces and the corresponding reduction of the potential barrier AEc = ^^^^^^^^^ 
will be much smaller then the effect induced by the channel A£^o ~ A£^i. 



5. Analytical approximation of the effect of the amplification 

A simple one-dimensional approximation of Poisson equation for a long (length d) and narrow 
(radius r) channel can be derived [U |2] and written on its dimensionless form as: 

= -/3A(1 - x) - aP(x) - a(5(x - xo) - a5(x - xi) (4) 
$(0) = 0, $(1) = A 

where B = — d?e_ -g ^ function of the effective dielectric parameter e = — j-ttj-tt- The 

factor a = K?d? where = y^fHfHa^^o^ ^^le Debye length with vq = nr'^d; the volume 
of the cylindrical channel. Ut = ksT is the thermal energy; where ks is the Boltzmann's 
constant and T the absolute room temperature. P{x) is the permanent charge distribution on 
the protein atoms; independent of the electric field. The 6 function represents two moving ions 
at different locations xq and xi on the channel axis. Because the Poisson equation is linear 
the contribution from the four terms in the right hand side can be considered independently. 
Here we concentrate on the interaction between two ions in the channel in the absence of any 
external field and therefore leave only the two last terms in the eq. (jH) (see, however, for the 
full analytical solution of the Poisson equation elsewhere). 



$(x) = - a G{x,s)P{s)ds 
Jo 



a 



sinh y/j3x sinh y/j3{xo — 1) + sinh y/j3x sinh y/j3{xi — 1), < x < xi 
, sinh y^xo sinh y^(2; — 1) + sinh sinh y^(a;i — 1), xi < x < xq 
V sinh V P [ sinh ^xo sinh ^{x - 1) + sinh ^xi sinh ^{x - 1), xq < x < 1 



(5) 



Where the Green function is given by: 



. _ 1 J sinh y^x sinh y^(s — 1), < x < s , . 

^"^'""^ ~ V(3smh,/(3[ sinh^ssinhV:5(x-l), s<x<l 

The distribution of the permanent charge is modelled using a narrow Gaussian distribution, 
consistent with the fact that the charge is mainly concentrated at the central part of the channel 
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6. Conclusion 

It was shown that the energy fluctuations due to the correlation of the motion between ions in 
strong electrolytes can be substantially enhanced in the open ionic channels with low dielectric 
constant. Therefore, ion channels can be taught of as electrostatic amplifiers of the charge 
fluctuations. This in turn may lead to an exponentially strong modulation of the potential barrier 
for the ion at the selectivity site. The resulting enhancement of the transition probabilities of 
ions through open ionic channels can be leading order contribution to the transition probabilities 
calculated using standard theories that neglect the effect of charge fluctuations. The effect of 
electrostatic amplification of the charge fluctuations may play substantial role in the nanoscale 
protein dynamics, including, in particular, protein folding. 
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